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Abstract. We provide an explicit formula for the Tornheim double series in 
terms of integrals involving the Hurwitz zeta function. We also study the 
limit when the parameters of the Tornheim sum become natural numbers, and 
show that in that case it can be expressed in terms of definite integrals of 
triple products of Bernoulli polynomials and the Bernoulli function Ak(q) := 

kC'(l-*,ff). 



1. Introduction 

The function 



1 

(1.1) T(a,b,c) = V V — — — — ■ -, a, b, c e C, 

^— ' ^— ' n a hi [n + m)° 

n=lm=l v 1 

was introduced by Tornheim in |25|. We provide here an analytic expression for 
T(a, b, c) in terms of the integrals 

(1.2) I(a,b,c)= f t(l-a,q)((l-b,q)ai-c,q)dq 

Jo 

and 

(1.3) J(a, b,c)= f C(l - a, q)((l - b, q)({l -c,l-q) dq. 



Here £(z, q) is the Hurwitz zeta function, 

(i.4) a^i) = E 



^ (n + q) z ' 

defined for z € C and q ^ 0, — 1, — 2, - ■ ■ . The series (II .4|) converges for Rez > 1 
and C(z, q) admits a meromorphic extension to the complex plane with a single pole 
at z = 1 as its only singularity. 

In the case where the parameters a, 6, c in l|l.l|) are positive integers, the Torn- 
heim sum can be expressed in terms of the Rlemann zeta function 



oo 

(1-5) CM = CM) = £- 



n=l 



its derivatives, and integrals related to the families (|1.2fl and (|1.3|) . as given in 
Theorem 11.11 below. 
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We use the notation 

(1.6) az,q) :=C(1 -*,<?), 

defined for z 7^ and q 7^ 0, — 1, —2, 

The results presented here are a continuation of I12| where we have provided 
many explicit evaluations of definite integrals containing £(z, q) in the integrand. 
For instance, if Re a > 0, then 



(1.7) 



C{a,q)dq = 0, 



and, for Re a > 1 , Re b > 1 , we have 



(1.8) 

and 

(1.9) 



C(a, q)C(b, q) dq = 2 [ 2 ( °™ C(« + 6) cos (|(a - 6)) 



C(a, ? )C(6, 1 - g) dq = 2 [ 2 ( °] a r + ( ? CO* + 6) cos (|(a + 6)) 



Lerch's evaluation [50] , 
(1.10) 



dz 



= lnT(g) -ln\/2^r, 



yields integrals involving the loggamma function. For instance, 
(1.11) 



Li = [ \nT{q)dq = lnV2^ 
Jo 



and 



(1.12) L 2 = f ln 2 T 
Jo 



(g)dg 



with 
(1.13) 



7 2 7T 2 7L1 4 

- — I h - — i + - 

12 48 3 3 



A = 7 + In 2n. 



Li 



IT 2 2tt 2 



We expect that the methods developed here will provide analytic expressions for 
the constant 



(1.14) 



L 3 = [ \n 3 T(q)dq. 
Jo 



The series T(a,b,c), for a, b, c e K — N, is given in Theorem 12.41 in terms of 
integrals ljl.2|l and (|1.3Jl . The evaluation of the Tornheim series for integer values 
of the parameters are expressed in terms of some definite integrals: 
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Theorem 1.1. The Tornheim sums T(nx, 112,113) can be expressed as a finite ex- 
pression of the Riemann zeta function, its derivatives and the integrals 

K m , n = [ ^- m Hq)B n (q)lnT(q)dq, 



K* = I ^- m \l~q)B n (q)\nT(q)dq, 



Z m ,u = / ^- m) (q)^- n Hq)\nT(q)dq, 
Jo 

Z* = [ ^- m \q)^- n \l-q)lnT(q)dq. 



Here B n (q) is the Bernoulli polynomial, T(q) is the classical gamma function and 

?/>^ m '(<z) = ip(—m,q), where 



(1.15) ^(z,q) = e-i*£- 



ryz CO + 1) g) 

e ' 



T(-z) 

is the generalization of the polygamma function introduced by the authors in |13| , 
with z an arbitrary complex variable. Some properties of ip{z, q) are given in 
Appendix 2. The closed form evaluation of the integrals in Theorem 11.11 will be 
discussed in a future paper. The proof of Theorem ll.ll is given in Sections and 21 

The series (11.1(1 converges for Re a, Re b, Rec > 1. Matsumoto showed that 
it can be continued as a meromorphic function to C 3 , with all its singularities 
located on the subsets of C 3 defined by one of the equations 

a + c = 1 - /, 6 + c = 1 - /, a + b + c = 2 with I E N . 

The literature contains many techniques to evaluate some particular cases of 
T(a, b 7 c). For instance, the case c = is evaluated simply as 

T(a,b,0) = C(a)C(&)- 

The elementary identity 

(1.16) T(a,b-l,c+l) + T{a-l,b,c+l) = T(a,b,c) 
and the symmetry rule 

(1.17) T(a,b,c) = T(b,a,c) 

has been used by Huard et al JTj to give the explicit expression 

i=l ^ a 1 ' i=l \ * ' 

in the case that both a and b are positive integers. Here N = a + b + c. If N is an 
odd positive integer greater than 1, then the sum T(i, 0, N — i) is evaluated as 



T(i,0,N-i) = {-!)* £ ( i-\ )CmC(N-2j) 



3=0 

\i/2\ 



3=0 
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The evaluation of T(i, 0, N — i) in the case N even remains open. The techniques 
introduced in this paper have allowed us to evaluate the sum T(a, 0, c) in terms of 
integrals similar to the ones discussed here. Details will appear in |14|. 
The multiple zeta value, also called Euler sums, are defined by 

(1-18) C(Ma,...,ifc) - J2^r4 — ~ k 

raj 1 n 2 ■ ■ ■ 

where the sum extends over m > n% > ■ ■ ■ > > 0. The sum T(a, 0; c) is £(c, a). 
A general introduction to these sums is provided in Chapter 3 of Borwein-Bailey- 
Girgenson 0. 

The identities of Tornheim |25| for T(n%, 712,713) are based on an elementary 
identity for series: let / be monotone decreasing and f{x) — > c as x — > 00 and 
define 

, n f(m) f(n) f{n + m) 
<P{n,m;f) = — — ■ r + — - — ■ . 

n(n + m) m(n + m) nm 

Then 

°° °° f(r) - c 

n,m— 1 r— 1 

The special case f(x) — 1/x yields 

(1.19) T(l,l,l) = 2C(3) 

and f(x) = \jx a ~ 2 produces the relation 

2T(a - 2, 1, 1) - T(l, 1, a - 2) = 2((a). 
Among the many evaluation presented in |25) we mention 

a-2 

T(l, 1, a - 2) = (a - l)C(a) - £ C(0C(« - *), 

1=2 

T(a- 2, 1,1) = ^T(l,l,a-2)+C(o) 



and 



T(l,0,a-l) = ir(l,l,a-2). 



Subbarao and Sitaramachandrarao 23 give 

(1.20) T(2n, 2n, 2n) = | £ f^" !_* 7 ^ C(^)C(6n - 2i), 

i=0 ^ ' 

and this is complemented by Huard et al. with 

(1.21) T(2n+ l,2n+ l,2n+ 1) = -4^ + 1 VO')C(6n - 2 * + 3 )- 

Boyadzhiev 0IH] has given elementary proofs of an expression for T(a,b,c) in 
terms of the function 



S(r,p) = £ 



OO rr{r) 
fin 



nP 

n=l 
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Here = l~ r + 2" r + . . . + n r is the generalized harmonic number. In [Hj the 
author establishes recurrences for the sums S(r,p) that permit to express them as 
products of zeta values in the case N = r + p odd. 

Tornheim double sums and other related ones appear as special cases of the zeta 
function Cg( s ) °f a semi-simple Lie algebra defined as 

C g (s) = ]Tdim(pr 
p 

where the sum is over all the finite-dimensional representations of g. Zagier |22I 
states that the special case g = s[(3) yields l|1.2U[) . The nomenclature for T(a, b, c) 
is not standard: Zagier (221 an d also Crandall and Buhler ^U] call T(a,b,c) the 
Witten zeta function. Tsumura |27| has evaluated some special cases of the sum 

oo 

W(p, q, r,s) = 



L^i mPni(m + n) r (m + 2n) s 

,71 — 1 ' 

under the parity restriction p + q + r + s is odd. This is the Witten sum corre- 
sponding to 50(5). 

These sums also have appeared in connection with knots and Feynman diagrams, 
see JHl I0r details. 

The Bernoulli function 

(1.22) A k (q) = kC(l-k,q), keN 

introduced in [T3|, plays an important role in the evaluations presented here. Adam- 
chik |2] proved the identity 

fc-i 

(1.23) C'(l-M) - ?(l-k)+Y,{-l) k - 1 - i 3lQj,k-i{q) lnr j+1 (g) 

j=o 

where 

(i-24) Q k Ai) = X>-«) n-J '(^) {i} 

is the Stirling polynomial and the generalized gamma function T n (q) is defined 
inductively via 

(1.25) r n+1 ( q + i) = ^±iM 

Ti(g) = T(q) 

r„(i) = l. 
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1>I OLcxLlOIl. 






r(«) 


is the gamma function, 






is the Riemann zeta function defined in (|1.5I). 




CM 


is the Hurwitz zeta function defined in 1)1.411. 




C(z,q) 


is a shorthand for £(1 — z, q), 




C±(z,q) 


denotes the combination £(z, q) ± £(z, 1 — q), 




ifj(z, q) 


is the generalized polygamma function, defined in 


Jl 1 r.l 

l|l.lo|) , 


^- n) (q) 


is the balanced negapolygamma function, defined 


in 


A n (q) 


is the Bernoulli function, defined in l|3.6[). 




B n {q) 


is the Bernoulli polynomial of degree n, defined in 


1731 and |Ol 


B n 


is the n-th Bernoulli number, 




H n 


is the n-th harmonic number, H n = 1 + ^ + • • • + 


l 

n ' 


h„ 


is a shorthand for H n -\, 




1 


is Euler's constant, 




A 


equals 7 + In 2ir, 




A± 


equals A 2 ± 





2. The main identity 

We now provide an analytic expression for the Tornheim double series T(a, b, c) 
in terms of the integrals l|1.2l) and l|1.3|) . The analysis of its behavior as the pa- 
rameters become integers is described in Section OH The proof employs the Fourier 
representation for £(z,q): 



(2.1) £(z,q) 



2I» 
(2tt) z 



cos ■ 
V 2 



/7rz\ cos(27rqn) 



£ 



7rz\ sin(27T(7n) 



valid for Re z > 1 and < q < 1 , given in 0] . 
From 1)2. 1|) we obtain 



(2.2) 

and 
(2.3) 

where 
(2.4) 



n=l 



2/.(*)£ 



cos 27rqn 



sin 27rqn 



n=l 



2r(z) /7TZ\ 

/cU) = -; — r— cos — and f s (Z 

Jc\ 1 y 2 J JS\ 



((z,q) + C(z,l-q) 

C(z,q)-az,l-q), 
2T(z) 



(27T) 



■ sin 



(t) 



For a function h(a, b, c) we denote 

(2.5) h sym (a, b, c) = h(a, b, c) + h(b, c, a) + h(c, a, b) 

and 

(2.6) h nsym (a, b, c) = -h(a, b, c) + h(b, c, a) + h(c, a, b). 
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Proposition 2.1. Let a, b, c G M. Then T(a,b,c) satisfies the relations 

(2.7) / c (a)/ c (&)/ c ( C )T^ m (a, 6, c) = /(a, 6, c) + J^ m (a, 6, c) 

and 

(2.8) f s (a)f s (b)f c (c)T ns y m (a,b,c) = I(a,b,c) - J ns ^ m {a,b,c), 
where f c and f s are defined in 1)2. 4[l . 

Proof. Multiply three series of cosine type in 1)2.2)1 to obtain that 

oo oc oo _^ „i 

8/c(a)/ c (&)/c(c) V V V — — t— I cos(2irqn) cos(2?rgm) cos(2?rgr) dq 
£ — ' * — ' * — ' 77 a 771 r c / n 

n— 1 m— 1 r— 1 

equals 

[C(<z, g) + C(o, 1 - g)] [C(b, q) + C(b, 1 - g)] [C(c, g) + C(c, 1 - g)] dg- 
The identities 

4cos(7t) cos(?j) cos(t«) = cos(?i + v + w) + cos(u + v — w) 

+ cos(t* — v + w) + cos(it — v — w) 

and 

(oifj^O 



U 



1 

cos(2-7rgj) dq 



1 if J = 0, 



reduce the left-hand side to 2f c (a)f c (b)f c (c)T sym (a, b, c). To complete the proof of 
the first identity, we expand the products of zeta functions to write the integral as 
a sum of eight different integrals, which can be reduced to the right-hand side of 
()2.7)l , by selectively performing the the change of variable q — > 1 — q in half of them. 

The second identity is obtained by considering the only other nonvanishing triple 
product integral, namely, that of sin(27rgn) sin(27T(7m) cos(27rqr). □ 

The case of Proposition l2~T1 in which the parameters are integers will be our main 
interest in this paper. When the argument z of the function £(z, q) is a positive 
integer n, this function reduces to a Bernoulli polynomial, 

(2.9) ((n,q) = -±B n (q). 

In this case, due to the reflection property of the Bernoulli polynomials, 

(2.10) B k (l -q) = (-l) k B k (q), 

the function £(77, 1 — q) reduces simply to £(n, q) up to a sign: 

(2.11) C(n,l-g) = (-l) n C(n,g), 

so that the J-type integrals reduce to /-type integrals. 

Unfortunately, since the functions f c (n) and f s (n) vanish for 77 odd and n even 
respectively, the identities of Proposition 12.11 for integer parameters (a, b, c) = 
(771, 772, 713) are trivial except only in two cases: 



8 



OLIVIER ESPINOSA AND VICTOR H. MOLL 



(1) m, 712,113 are all even, 

(2) ni,n2 are both odd, and 713 is even. 

The first case is of special interest. It appears in as the reciprocity relation for 
a class of Tornheim series. 

Corollary 2.2. Let n±, ri2, 713 G N be even. Then 
(2.12) 

(n 7r \ni+n2+n3 

{n 1> n 2> n 3 ) ( L) 2(m - l)!(n 2 - l)!(n 3 - 1)! 1 X ' 2 ' 3j ' 

Corollary 2.3. Let n e N. Then 

(2.3, T(2„,2,2„, . i(-l)»( 2 ,r E( 4 " 2 ;l\- (2^-t t) r 



fc=0 

Proof. Corollary EH and (23 yield 

T(2n,2n,2n) = (-l)" +1 ^-i^ jf B 2n ( ) 3 dq. 

The value of the integral is given by Carlitz It can also be obtained directly 
from the formula for B n (q) 3 given in Appendix 1. □ 

Formula H2.13|) agrees with formula H1.20JI on account of the relation 

(2 - 14) C(2fc) = p ' 

valid for k e N . 

We now present an analytic expression for the Tornheim double series, valid for 
non-integer values of the parameters. 

Theorem 2.4. Let a,t,c€l and define 

r(i-z) tt 



(2.15) X(z) 

y ' w (27T) 1 " 2 (2tt) 1 - 2: F(z) siiiTrz 

For a, 6, c ^ N we have 

(2.16) T{a, 6, c) = 4A(a)A(6)A(c) sin(7rc/2) cos (f (a - 6)) [J(c, a, 6) + J(c, 6, a)] 

- cos (f (a + 6)) [I(a, b. c) + J(a, b, c)\ 
Proof. The difference of the two expressions stated in Theorem 12 . II yield 
(2.7, 2 Ucma , „, c) . (-^J^ - jj^jj) U(a. >,c ) + J(a, „, c)] 



1 1 



[J(c, a, 6) + J(c, 6, a)] 



Jo(o)/c(6) fs(a)f s (b), 

and the result follows directly from here. The values of a, 6, c € N are excluded due 
to the singularity of X(z) for zeR □ 
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3. The limiting case 

The goal of this section is to analyze the result of Theorem 12. 41 as the parameters 
a, b, c approach positive integer values. The notation a = n± + £i, b = n% + £2, c — 
ri3 + £3 with rij £ N and Ej — ^ is used. 

We start by writing 

(2ir) a+b+c 

with, according to l|2.4[) and l|2.f 71) , 



(3.2) f(a,b,c):=- 

Cc 



J(a, b, c) + J(a, 6, c) + J(c, a, 6) + J(c, 6, a) 



C a Cb 

-/(a, 6, c) — J(a, b, c) + J(c, a, b) + J(c, 6, a) 



S a Sb 

where 

c a = cos(7ra/2) and s a = sin(7ra/2) 
and Ch, Sf,, c c , s c are similarly defined. 

Our first task will be to obtain the limit of c c T(a, b, c) as both a and b approach 
positive integer values. The functions c a and s a have the property that, as the 
real number a approaches an integer, one of them tends to zero while the other 
tends to +1 or —1, depending on the parity. Explicitly for n an integer and e an 
infinitesimal quantity, 

{(— l) n / 2 + o(e), n even 
(-1) 2 _ £ + ( £ ) ) n odd, 

J(-l)" /2 5e + o(£), neven 

(3-4) S n +e = < 

[(—1) 2 +o(e), n odd, 

Hence, in order to compute the limit we are seeking, we need to expand the nu- 
merators inside square brackets in H3.2JI up to order £i£2- This is accomplished 
by replacing both £(a, q) and £(6, q) by their Taylor series expansions around an 
integer value of its first argument, 

C(n + e,q) = £{n,q) +e('(n,q) + o(e) 

(3.5) = --[B n (q)+eA n (q)] + o(s), 

n 

according to l|2.9|l and the definition of the Bernoulli function 

(3.6) A k (q):=k('(l-k,q), 

studied in ^J^]. For instance, using the shorthand notation 

(3.7) </(<?)}:= / f(q)dq 

Jo 
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we have 

(3.8) J(c,a,b) 



a = m + £1 
b = ri2 + £2 



nin 2 



[(B ni {q)B n2 {l-q)C{c,q)) 



+ e 2 <-Bni(«)A 2 (l-«)C(c.g)) 

+ £i£2 (A ni (q)A n2 (1 - g)C (c, g))] + o(£i£ 2 ). 

Using the reflection property 1|2.1U|) of the Bernoulli polynomials and the invariance 
of the integration (|3.7|l under the change of variable q — + 1 — q, we find the following 
result for the numerators inside square brackets in (|3.2[1 (the upper sign corresponds 
to the numerator of the first term and the lower sign corresponds to the numerator 
of the second term): 

(3.9) ± I(a, b, c) ± J(a, b, c) + J(c, a, b) + J(c, b, a) = 

— { ± [±1 + (-1)" 1 ] [±1 + (-1)" 2 ] (B m (q)BnMC(c,q)) 

+ ei [±1 + (-1)" 2 ] (A ni (q)B n2 (q)(+ (c, g)> 
+ e 2 [±1 + (-1)" 1 ] (B ni (g)A„ 2 (g)C+ (c, g)) 

+ ei£ 2 [± (A ni (q)A n2 (g)C+ (c, g)) + < A ni (g) A n . 2 (1 - g)C+ (c, g)>] } 

+ o(eie 2 ). 

where 

(3.10) c+(«,?):=C(«,9)+C(2,i-g). 

We now examine the behavior of the Tornheim sum as £i,£ 2 — > 0. The limiting 
value is obtained from (|3.3|l . I|3.4|) and l|3.9|) . Observe that T(a,b,c) = T(b,a,c) so 
only three cases are presented. 

Theorem 3.1. Suppose m, «2 £ N and c € K\N. Then the Tornheim double 
series T{n\, n 2 , c) are given by 



v ' ; v 7 4m!n 2 ! T(c) cos(ttc/2) 



B ni (q)B n2 (q)((c,q) dq 



1 



-a / ^ni(g)^n 2 (g)C+ (c,g)cig + — / A ni (q)A n2 (1 - g)C+ (c, g) <2g 



for rii, n 2 even; 



(3.12) (-1) s V y 



(2tt) c 



4ni!n 2 ! r(c) cos(ttc/2) 



^ / L'„, U/U„,U/k-. {C.q)<lq 

+ - I A ni (q)B n2 (g)C+ (c, g) rfg 



for ni even and n 2 odd, and 
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(3.i3) 



4rii!n 2 ! T(c) cos(7rc/2) 



/ B ni (q)B n2 (q)((c,q)dq 
Jo 



~ ~2 I An i (l) A n2 («K+ (c, q) dq - \ ( A ni (q)A n2 (l-q)( + (c,q)dq 
Jo ™ Jo 

for m, n 2 odd; 

The final step in the process is to let c = 713 + £3 and let £3^0. For n even we 
simply have 

-^ T =--(-D" /a B B («), 
cos(7rn/2) n 

whereas for n odd, 

lim it/ C 'l = -^(-1)^-^(9) + A,(i - 5)] • 

c^n C0S(7TC/ 2) 71 TT 

The value of T(ni, n 2 , 7J3) is thus expressed in terms of integrals of triple prod- 
ucts of the Bernoulli polynomials Bk{q) = —k((l — fc, g) and the function Afc(g) = 
fcC'(l-fc,g). 

Define the following families of integrals: 

(3.14) Ri(ni,n 2 ,n 3 ) = B ni (q)B n2 (q)B n3 (q) dq 

Jo 

1 f 1 

(3.15) R 2 {ni,n 2 ,n 3 ) = - / B ni (q)B n2 (q)A n3 (q) dq 

n Jo 



(3.16) R3(ni,n 2 ,ns) = — I A ni (q)A n2 (q)B n3 (q) dq 







1 r 1 

(3.17) R 4 ( ni ,n 2 ,n 3 ) = — A ni (q)A n2 (l - q)B n3 (q) 

n Jo 

1 f 1 

(3.18) i? 5 (ni,n 2 ,n 3 ) = — / A ni (g)A„ 2 (g)A„ 3 (g) dg 

Jo 

1 r 1 

(3.19) R 6 (n 1 ,n 2 ,n3) = —» / A ni (g)A„ 2 (g)A„ 3 (1 - g) dg. 



These integrals are all symmetric under interchange of their first two arguments 

(ni and n 2 ), except for R4 which is antisymmetric if n 3 is odd. 

Define 



(3.20) p(n) 



'(-1)" /2 , 



n even 
(— l) 11 ^ - , n odd. 



Theorem 3.2. Let a = ni+n 2 +n 3 . Then the Tornheim double series T(m,n 2 ,n 3 ) 
is given by 

(3.21) T{ ni ,n 2 ,n 3 )=p(a) - ^\ — 7 T R (n u n 2 , n 3 ) 

2ni! n 2 ! n 3 ! 
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where Tn(n\, n 2 , n 3 ) can be expressed in terms of the functions Rj : 1 < j < 6 as 
follows: 



Case 1. ni, n 2 and n 3 are even: 

(3.22) r R (n 1 ,n 2 ,n 3 ) = --#i(ni,n 2 ,n 3 ) + R 3 (n 1 ,n 2 ,n 3 ) - R 4 (n 1 ,n 2 ,n 3 ) 

Case 2. n\ and n 2 are even; 713 is odd: 

(3.23) T R (ni 1 n 2l n 3 ) = -# 2 (ni, n 2 , n 3 ) + #5(711, rc 2 , n 3 ) 

+ R 6 (ni,n 2 ,n 3 ) - Re(n 3 , ni, n 2 ) - R e (n 3 ,n 2 ,ni) 

Case 3. ni is even, n 2 is odd, and 713 is even: 

(3.24) T R (m,n 2 ,n 3 ) = -# 2 (n 3 ,m,n 2 ) - i? 2 (n 3 , n 2 , n x ) 
Case 4. ni is even; n 2 and 713 are odd: 

(3.25) Tfl(m, 712,713) = R 3 {n 3 ,ni,n 2 ) + # 3 (n 3 , n 2 , ni) 

+ #4(711, n 3 ,7J 2 ) + -R 4 ("-2,«3,«l) 

Case 5. ni and n 2 are odd; n 3 is even: 

(3.26) Th(tii, 712,713) = --#1(711,712,713) + # 3 (ni,n 2 ,n 3 ) + #4(711, n 2 ,n 3 ) 
Case 6. m, n 2 and 713 are odd: 

(3.27) T R (ni,n 2 ,n 3 ) = -R 2 (ni,n 2 ,n 3 ) +#5(711,712,713) 

+ #6(ni, 712, 713) + #6(713, 711,713) + #6(713, n 2 ,m). 

The closed form evaluation of the Tornheim sums T(ni,n 2 ,n 3 ) has thus been 
reduced to that of the integrals Rj . A partial evaluation of these integrals is pre- 
sented in the next section, in terms of new family of integrals, closely related to 
Rj. 

4. A NEW FAMILY OF INTEGRALS 

The evaluation of the integrals Rj is most conveniently organized in terms of a 
new family of integrals Qj , defined in terms of the balanced generalized polygamma 
function, introduced in |13j . by 

1 

B ni (q)B n2 (q)B n3 (q) dq 



(4.1) 


Qi(ni. 


n 2 . 


n 3 ) 


(4.2) 


£2(711, 


n 2 . 


n 3 ) 


(4.3) 


Q3(ni : 


n 2 . 


n 3 ) 


(4.4) 


Qi{n 1: 


n 2 . 


n 3 ) 



(I 



B ni (q)BnM^ n3 \q)dq 











1 



i>(-^\q)^- n *\l-q)B n3 (q)dq 
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(4.5) Q 5 (n 1; n 2 ,n 3 )= / (q)^- n *\q)i,^ (q) dq 

Jo 

(4.6) Q 6 (n 1 ,n 2 ,n 3 )= [ {q)^- n -\q)^-^\\ - q) dq. 



These integrals satisfy the same symmetry properties as their i?-analogs. In ad- 
dition, the Q-integrals satisfy homogeneous recursion relations, which allow their 
evaluation in terms of a few basic ones. The relation among the families Rj and 
Qj using the identity 

(4.7) A m (q) = m\^- m \q) + h m B m (q) 

is given in Appendix 3. 

In this section we present recurrences for the integrals Qj . The initial conditions 
require a variety of definite integrals listed below: 



^- m \q)B n {q)dq, 
^ m \q)^- n \q)dq, 



• The integrals 




(4.8) 




(4.9) 


M ro ,„ = 


and 




(4.10) 


M* = 

m,n 


which will be evaluated in 


• The families of integrals 


(4.11) 


K — 


(4.12) 


K* = 

m , n 


(4.13) 


Zm.n — 


and 




(4.14) 


ry* 

m,n 



1 

lp { - m) {q)lp ( - n) {l -q)dq, 



^- m \q)B n (q) Inr(g)dg, 
^—)(1 - q)B n {q) \nT(q)dq, 

l 

^- m \q)^- n \q) \nF{q)dq, 



V> ( ~ m) (q) i> { ~ n) (1 - q) In Y(q) dq. 
The closed form evaluation of these functions is left as an open question. 
The integral Q\. The explicit value of Qi(ni, 112,113) was given by Carlitz 9 : 

L(m+n 3 -l)/2j r / \ / n 

(4.15) Qi(n u n 2 ,n 3 ) = (-l)" 3+1 n 3 ! V ^1 J + n 2 W 

(ni+n 2 -2fc-l)! 
X (n 1 +n 2 +n 3 -2A ; )!^ 2fe ^ 1+ " 2+n3 - 2 ' £ - 
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The integral Q 2 . This integral is obtained directly from the formula given in 



(4.16) 



((z + 1, q)^{z', q)dq = 2(2tt)*+* 'r(-z) { ~C(-3 - z') sin '-{z - z') 



( 7 + In 2tt)C(-z - z') - C'(~« ~ z')J cos -(z - 



valid for Rez, Rez' < and Re(z + z') < — 1. The evaluation at z 
z' = —n, with to, n G N gives 
-l 



-to and 



(4.17) / B m {q)^- n \q)dq 



2m! 



| — £ (to + n) sin — (to — n) 



(27r) m+n 

7T ~1 

— [(7 + In 2n)( (to + n) — £ (to + n)] cos — (m — n) | . 

The evaluation of Q 2 (ni, n 2 , TI3) follows from 14.17|l and the representation (|7.3|1 
for the product of two Bernoulli polynomials: 

Theorem 4.1. Let ni,n 2 , n$ G N and let a = n-i + ?^ 2 + "3- Then 
(4.18) 



?2(ni,n 2 ,n 3 ) = 2(-l) n3 ^ 



fc=0 



»i 



n 2 



(m + n 2 -2fc-l)! 



(27T) 



a-2k 



x {-l) k B 2k {-sin — C (a - 2k) - cos — [(7 + m2^)C (a - 2k) -('(a- 2k)]} 

where k(ni,n 2 ) = Max{ |_ii/2j , |_n.2 /2J } . The constant term in (|7.3I) gives no con- 
tribution on account that the function ^~ n ^(q) is balanced for n G N. The integrals 
Q 2 are only needed in the case that a = n\ + n 2 + n% is odd, equal to 2N + 1, say. 
In this case, sin(7ra/2) = (— 1)^ and cos("7ra/2) = 0, so that 



(4.19) Q 2 (ni,n 2 ,n 3 )=7r(-l) Ar+ " 3 

k(ni ,712) 



E 

k=0 



n 2 



(nx + n 2 -2k-l)\, k 



(27T) 



i-2k 



{-l) k B 2k C{a-2k). 



Proof. The details are elementary. 

The integral Q3. We now produce a recurrence for 



□ 



Qa(ni,n 2 ,n 3 ) = 

The basic tools are the relations 

d 



^-^(q)^-^(q)B n3 (q)dq. 



(4.20) 
and 
(4.21) 



dq 



^S m (q) = mB ro _i(g), 



valid for m G No, and the fact that both the negapolygamma functions and the 
Bernoulli polynomials are balanced for the range of indices we wish to consider, 
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i.e., ^- m \l) = </> ( ~ m) (0), for m > 2, and B m (l) = B m (0), for all m. 

Theorem 4.2. Let n\, n 2 , n 3 £ N with n\,n 2 > 1. Then 

(4.22) (n 3 + l)<2 3 (ni, n 2 , n 3 ) = -Q 3 (ni - 1, n 2 , n 3 + 1) - Q 3 (ni, n 2 - l,n 3 + 1). 
Proof. Start with 

(n 3 + l)Q 8 (ni,n2,n3) = j\ { - ni) (q)^- n2] (q)-^B n3+1 (q)dq, 
integrate by parts and observe that there is no contribution from the boundary. □ 

The recurrence shows that the value of Q 3 (fT-i, n 2 , n 3 ) can be obtained from the 
values of 

(4.23) Q 3 (l,m,n) = K m , n + C'(0)JV m>n , 

in view of t/* (?) = m r(g) + C'(0)j the symmetry of the integral Q 3 under in- 
terchange of its first two arguments, and the definitions l|4.11|l and (|4.8|l of the 
integrals K m . n ,N mtn . 

The integral Q 4 . Similarly, for ni, n 2 > 1, we have that 

Q 4 (n 1 ,n 2 ,n 3 ) - f\/- n ^ (q)^~ n ^ (1 - q)B n3 (q) dq 



J T13 \ 
JO 

satisfies the recurrence 

(4.24) (n 3 + l)Q4(ni,n 2 ,n 3 ) = -Q4(ni - l,n 2 , n 3 + 1) + Q 4 (ni,n 2 - l,n 3 + 1), 
so that it can be obtained from 

(4-25) Q 4 (l,m,n) = i^ n + C'(0)(-l) n iV ro , n , 

where is defined in (|4.12|) . 

The integrals Q5 and Q6- Similar arguments show that for n\, n 2 > 1 we have 

Q6(ni,na,ns) = f ^-^\q)^- n *\q)^-^) ( q ) dq 
Jo 

Q 6 (ni,n 2 ,n 3 ) = f ^~ n ^ (q)^~ n ^ (q)^~ n ^ (1 - g) 



./() 

satisfy the recurrences 

(4.26) Qs(m,n 2 ,n 3 ) = -Q 5 («i - l,n 2 ,n 3 + 1) - Q 5 («-i,"2 - l,n 3 + 1) 
and 

(4.27) Q6(ni,n 2 , n 3 ) = Q 6 (ni - 1, n 2) n 3 + 1) + Q 6 (ni, n 2 - 1, n 3 + 1). 

The initial conditions are 

(4.28) Q 5 (l,m,n) = Z m ,„ + C(0)Af m , n , 

(4.29) Q 6 (l,m,n) = Z^ n + C(0)M^ >n , 

where Z m , n , Z* l n , M m , n and M* l n are defined in (|4~T^1> . gHU), and (f^TTUf) . 

respectively. 
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5. Evaluation of integrals 

5.1. The product of a Bernoulli polynomial and a balanced negapolygam- 

ma. The explicit evaluation of the integral N mjl defined by l|4.8(l was obtained in 
(|4.17|l as a byproduct of the evaluation of Qi\ 

2n\ r 7r 7r 
(5.1) N m!n = , 2 s m+n \ 2 C (m + n) sin -(n - m) 

— [(7 + ln27r)£ (m + n) — C ( m + n )] cos — (n — m) > . 



5.2. The product of two balanced negapolygammas. The integral of the 
product of two balanced negapolygamma functions has been given in |12| : 
Let k, k' e N, k+ = k + k' and fc_ = k - k'. Then 

(5.2) M Kk ,= \ ^- k \q)^- k '\q)dq 

Jo 

2 



(27r) fc + 



[A + ak + )-2AC(k + )+C"(k + ) 



where A = 7 + In 2ir and A± = A 2 ± 7r 2 /4. 



Similarly, one finds 



M. 



k.k' 



(5.3) 



^- k) {q)^- k,) {l-q) dq 
(2^{ C ° S (I K ) [ A - C ^ - 2A C'( fc +) + ^"( fc +) 



+ 7rsin(|fc + ) [AC(fc + )-C'(fc+)]}- 



We have been unable to obtain closed form results for the integrals K m n n , 
Z m ,n and Z^ n , which involve the kernel lnT^) and one or two negapolygamma 
functions. 



6. Some examples 

In this section we describe some explicit evaluations of Tornheim sums. It is 
convenient to introduce the function 

(6.1) U m , n = f ^~ m \q)B n {q) lnsin(Trg) dq. 

JO 

The identity 

(6.2) In T(q) + lnT(l - q) = In tt - In sin(Trg) 
yields the relation 

(6.3) K m , n + {-l) n K* m ^ n = lnnN m , n -U m . n . 
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Example 6.1. We consider the value of T(l, 1, 2). This corresponds to Case 5 of 
Theorem El It yields 

1, 2) + .R 3 (1, 1,2) + i? 4 (l,l,2) , 



(6.4) 



T(1 ! 1,2) = 4tt 4 



and in terms of the Q 3 ; -family, 



(6.5) 



T(l,l,2) = 4rr 4 



The Qj -integrals are given by 



-g 1 (l,l,2) + ^Q 3 (l,l,2) + ^Q 4 (l,l,2) 



1 



Qi(l,l,2) = 

180 

Q 3 (M,2) = K 1>2 + C'(0)iVi,2 

Q 4 (l,l,2) = ^,2 + ^(0)^1,2 

and using the values C'(0) = — \ ln(27r) and N\p — C(3)/47r 2 we obtain 

T(l,l,2) - 4n 2 (K\ 2 + ^r !2 )-C(3)ln(27r)-l7r 4 . 
In terms of the [/-function this can be written as 



(6.6) 



7(1,1,2) 



-4ir 2 U h2 - C(3) In 2 



90 



-IT 4 . 



The identities ifTToj) and ((TT7|) yield T(l, 1, 2) = 2T(0, 1,3), and the method of 
Huard et al. |17) yields the values of T(ni,n 2 ,n 3 ) for N — ni + n 2 + n 3 odd and 
also for N — 4 and 6. For instance T(0, 1,3) = -|C(4)- This yields an evaluation 
of an integral of type U mjn : the value T(l, 1,2) = £(4)/2, the identity ^~ x \q) = 
lnr(g) +C'(0) and 



/ B 2 {q) ln(sin7rg) dq = - 
Jo 



03) 
2tt 2 



given in Example 5.2 of ^5 produce 



(6.7) 
It follows that 
(6.8) 



B 2 (q) In T(q) hi(sm.irq) dq 



7T 

240 



ln(47r)C(3) 
4tt 2 



7T 

240 



ln2C(3) 
4tt 2 ■ 



Example 6.2. The explicit expression for R 2 (ni, n 2 , 77,3) that can be obtained from 
(19.21) through l|4.15|l and (|4.18() permits the evaluation of T(ni, n 2l n-i) in the case 
ni, ^3 even and n 2 odd. For example 

T(2,l,2) = yC(3)-|c(5), 

T(2,3,2) = -^C(5) + 2C(7), 
6 

and 

T(4,3,2) = ^C(5) + yC(7)-^C(9). 
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Example 6.3. Define w = a + b + c to be the weight of the sum T(a, b, c). The 
results of the procedure described above for sums of small weight are given below. 



weight 3: 

m, <*, + nz U - C(3) + 1„2. (* 4 - «<M + 

weight 4: 

T(l,l,2) = -47r 2 C/ 1 , 2 -^-C(3)ln2, 
T(l,2,l) = -47r 2 (£/ 1;2 + 2[/ 2; i) , 
weight 5: 

T(l, 1, 3) = -8tt 2 + 2Z h3 + 2Z; i3 + AZ*^) - ( (5) 

/tt 2 A 2 tt 2 A ^ AC (A) A 3C ; (4) 2(" (4) 

T(1,2,2) = ^C(3)-|C(5), 

T(l, 3, 1) = -8^ 2 (^3 + 2Z 3 ,i + 2Z* 3 + 4^0 + y C (3) - 2( (5) 

ln2 '^_^_^_4C^)A 3H4) + 2C(4) 

1 45 30 90 7T 2 7T 2 7T 2 



2 



T(2, 2, 1) = 32tt 2 Z 2 , 2 + y C (3) - 3C (5) 

ln2W ^_^ + 4H4)A_2q4) 

45 180 7T 2 7T 2 

We have produced some partial results in the evaluation of the integrals K m>n , 
n and Z m ^ n , Z^ L n . These suggest that the value of the Tornheim sums can be 
expressed in terms of a small number of definite integrals. For instance, for m > 3 
odd, we have the relation 

n +i + In V2n(N m , n + m-/V m _i ;n+ i) 
^ B m (q)^(q)^- n -^(q) dq, 

which reduces the value of K m ^ n to that of Ki im+n _i plus the moments of the 
product ip{q)i)^~^(q)- Details will be presented elsewhere. 



7. Appendix 1: The Bernoulli polynomials 
The Bernoulli polynomials B n (q) defined by the generating function 
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The Bernoulli numbers B n = B n (0) satisfy 
(7.2) 



k=0 ^ ' 



For n > 1 we have the differential recursion B' n (q) = nB n _i(q) and the symmetry 
rule B n (l -q) = (-l) n B n (q). In particular, B„(l) = B n (0) for n > 1. 

The Bernoulli polynomials {B (q), Bi(q), ■ ■ ■ , B n (q)} form a basis for the space 
of polynomials of degree at most n. Thus the product B ni (q)B n2 (q) is a linear 
combination of Bj(q) for j = 0, • • • , ri\ + n^. It is a remarkable fact that this 
combination has the explicit form 



(7.3) B ni (q)BnM = E 



k=0 



n 2 



Bik 



-2k 



-B 



n 1 + n 2 -2k(q) 



, n n 1 + l "l'"2! „ 

+ l j (m+n 2 )! " 1+ " 2 ' 

where fc(ni, 712) = Max{ [«i/2j , [1-2/2] }. In terms of rescaled Bernoulli polynomials 
and numbers, defined by 



(7.4) 



B n (q)= B ^f-, S„ = S n (0) = ^, 
n! n! 



relation (|7.3|l has the simpler form 

k{n\ ,n 2 ) 

(7.5) ^(5)^(9)= J2 [(^Z-f' 1 ) 

k=0 



(n 1 +n 2 -2k-l\ 
{ na-1 1 



B2kB ni + rl2 



-2k{q) 



+ (-l)" 1 + i i?„ 1+ „ 



In theory, l|7.3[l yields expressions for a product of any number of Bernoulli 
polynomials. For example, 



L»/2J 



(7.6) 
or 

(7.7) 

and 
(7.8) 

Bl{q) 



Bl(q) = E 



n \2k) p r („\\( 1 \n+l -^ 2 " 
£>2k£>2n-2k(q) + [— J-J 72^- 



fc=0 



L«/2J 



( 2 ;) 



5^(9) 



2 E ( 2n ^i~ 1 )^2 fe S 2 „_ 2fe (q) + (-l)"+ 1 J B 2n , 



n—k 



L»/2J /n\ 

E E 

n — k '-^ 

k=Q j=a 



[2n - 2k 



+ (-1)' 



n+l 



^„( g) + 2n! 3 E 



+ 2(n-fc) 

L"/2J 



-°2j B3 n -2k-2j(q) 

in -2k- 2j 



2n-2k - 1\ B 2 kB 3n -2k 



n — 1 



J (2fc)!(3n-2fc)! 
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Ln/2j 



fe=0 



(7.9) Bl{q) = 2Y J ( 2n ;l\- 1 )B2 k 

n — k 



(3n-2k-2j-l\ 
\ 2n-2fc-l J 



[n/2\ 



fc=0 



3=0 

Integrating the relation (|7.1|l yields 

(7.10) [ B n {q)dq = 0, for n > 1 

Apostol |3j gives a direct proof of 



Bij B 3n _2k-2j(q) 



B 2n B n {q) + 2 (^n-t 1 ) 



B ni {q)BnMdq = ni!n2 ' 



, |"ni+H2i 



(7.11) , , , , .j , , • ' , ,1 

Jo (ni+n 2 y. 
for m, n 2 G N. 

The Bernoulli polynomials appear also as special values of the Hurwitz zeta 
function 



(7.12) 



8. Appendix 2: The generalized polygamma function i>{z,q) 
The polygamma junction is defined by 

(8.1) 

where 

(8.2) 



^)=^°)(g) = f hir(g) 

is the digamma function. 

The function tpt" 1 * 1 is analytic in the complex q-plane, except for poles (of order 
m+1) at all non-positive integers. Extensions of this function for m a negative inte- 
ger have been defined by several authors ^ ^1 El • These are the negapoly gamma 
functions. For example, Gosper defined 

V>_i(g)=lnr( g ), 



(8.3) 



^-k{q) 



iP- k+ i{t)dt, k>2, 



which were later reconsidered by Adamchik pQ in the form 
(8.4) 



i>-k(q) = TT^yy J (q- t) k ~ 2 hxT(t)dt, k>2. 



These extensions can be expressed in terms of the derivative (with respect to its 
first argument) of the Hurwitz zeta function at the negative integers ^El- The 
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definition (jHSJ 

can be modified by introducing arbitrary constants of integration 
at every step. This yields different extensions differing by polynomials: 

satisfying 

Pn(q) = -rPn+i{<i)- 
aq 

A new extension of ip^ m \q) has been introduced in |12|. in connection with inte- 
grals involving the polygamma and the loggamma functions. These are the balanced 
negapoly gamma functions, defined for mgNby 

(8.5) V ( " m) (<z) := —. [A m (q) - H m _iB m {q)] . 

ml 

Here H r = l + l/2 + -- -+l/ris the harmonic number (H = 0), B m (q) is the m-th 
Bernoulli polynomial, and 

(8.6) A m (q)=mC(l-m,q). 

A function f(q) is defined on (0, 1) is called balancedii its integral over (0, 1) vanishes 
and /(0) = /(l). In ^2] we have shown that 

(8.7) 4^ (-m) 0?) = ^ ( ~ m+1) (g), m e N. 
dq 

The function tp(z, q) defined in (|1.15l) represents an extension of these polygamma 
families to q £ C. Its main properties are presented in the next theorem. The de- 
tails appear in |13| . 

Theorem 8.1. The generalized polygamma function ip(z,q) satisfies: 

• For fixed q € C, the function ip{z, q) is an entire function of z. 

• For m e Z : ip(m, q) = ^ {m ~>{q). 



• It satisfies 



d 

—%p(z 7 q) = ijj(z + l,q). 



9. Appendix 3: The relation between Q 3 and Rj 

Using the relation (|4.7|l we can express the integrals Rj, defined by l|3.14|l - 
(|3.19|l in terms of Q j: defined by (|4.1() - H4.6(I . Recall that, for n > 1 we have 
h n = 1 + 1/2 H h l/(n - 1) and /ii = 0. 

(9.1) Ri(n 1 ,n 2 ,n 3 ) = Qi(n 1 ,n 2 ,n 3 ), 

(9.2) irR 2 (ni,n 2 ,n 3 ) = n 3 \Q 2 {n 1 ,n 2 ,n 3 ) + h n3 Qi(rii,n 2 ,n 3 ), 

(9.3) ir 2 R 3 (ni,n 2 ,n 3 ) = n 1 \n 2 \Q 3 {n 1 ,n 2l n 3 ) + n 2 \h ni Q 2 (n-i, n 3 , n 2 ) 

+ n 1 \h n2 Q 2 (n 2 ,n 3 , m) + h ni /i„ 2 Qi(ni, n 2 , n 3 ), 

(9.4) TT 2 R 4: (ni,n 2 ,n 3 ) = niln^.Q^nx,^,^) + (-l) n2 ni\h n2 Q 2 (n 2 , n 3 , nx) 

+ {-l) ni+n3 n 2 \h ni Q2(n 1 ,n 3 ,n 2 ) 
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(9.6) 7r 3 i? 6 (ni,n 2 ,n 3 ) 



(9.5) ir 3 R 5 (ni,n 2 ,n 3 ) 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



(-l) n2 h ni h n2 Qi(n l7 n 2 ,n 3 ) 7 

niln 2 \n 3 \Q 5 (n 1 ,n 2 , n 3 ) + iii\n 2 \h n3 Q 3 {n ll n 2l n 3 ) 
ni\n 3 \h n2 Q 3 (ni, n 3 , n 2 ) + n 2 \n 3 \h ni Q 3 (n 2l n 3l m) 
n 1 \h n2 h ri3 Q2{n2,ns, ni) + n 2 lh ni h n3 Q 2 (m,n 3 , n 2 ) 
n 3 \h ni h n2 Q 2 (ni , n 2 , n 3 ) + h ni h ri2 h n3 Q 1 (ni , n 2 , n 3 ) , 
ni!n 2 !n 3 !(5 6 (ni, n 2 , n 3 ) + (-l) n3 n 1 \n2\hn s Q3(ni,n 2 , n 3 ) 
ni!n 3 !/i„ 2 Q 4 (ni, n 3 , n 2 ) + n^.n^.h^Q^n-z,^, ni) 
(-l)" 3 ni!/i„ 2 /i„ 3 <5 2 (n2,n3,«i) 
(-ly^n^.h^h^Q^nx, n 3 , n 2 ) 



+ 



+ 




h n2 Q2(ni,n 2 ,n 3 ) 
Qi(ni,n 2 ,n 3 ). 
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